Hidden Order and Dynamics in Supersymmetric Valence Bond Solid States 
- Super-Matrix Product State Formalism - 
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Supersymmetric valence bond solid models are extensions of the VBS model, a paradigmatic model of 'solv- 
able' gapped quantum antiferromagnets, to the case with doped fermionic holes. In this paper, we present a 
detailed analysis of physical properties of the models. For systematic studies, a supersymmetric version of the 
matrix product formalism is developed. On ID chains, we exactly evaluate the hole-doping behavior of various 
physical quantities, such as the spin/charge excitation spectmm, superconducting order parameter. A general- 
ized hidden order is proposed, and the corresponding string non-local order parameter is also calculated. The 
behavior of the string order parameter is discussed in the light of the entanglement spectmm. 



I. INTRODUCTION 

Valence bond solid (VBS) models introduced by Affleck, 
Kennedy, Lieb and Tasaki''- are exactly solvable models that 
exemplify the gapped ground states in integer-S* spin chains 
conjectured by Haldana^ii. Though the VBS states, which 
are the exact ground states of the VBS models, are disor- 
dered spin liquids in the sense that their spin-spin correla- 
tions are exponentially dumped with a very short correlation 
length, there still exists a certain kind of "hidden order" cap- 
tured by the non-local string parameteri^. The existence of 
the hidden order highlights the exotic features of the Haldane- 
gapped antiferromagnets which are considered as manifesta- 
tion of the topological order of quantum spin chain^^. With 
recent increasing interests in the topological states of matter 
spurred by the discovery of topological insulators [See Reffisl 
for instance as a review], the VBS model and its variants are 
attracting renewed attention. Since the VBS states enable 
us to calculate many interesting quantities exactly, they of- 
fer a rare theoretical playground for the study of topological 
states of matter Due to their peculiar features, the VBS-type 
states have been investigated in a wide variety of contexts 
like quantum informatio 



topological orderi^iii, entan- 
glement entropyi^—, higher symmetric generalizations^^, 
and topological phase transitions22i^. 

In this paper, we present a detailed analysis of the re- 
cently proposedii supersymmetric generalization of valence 
bond solid (sVBS) states, The sVBS states are a precise 
mathematical realization of Anderson's scenario of high-Tc 
conductivity^^- and the idea of symmetry unification of super- 
conductivity and antiferromagnetism^l The sVBS states are 
hole-pair doped VBS states containing both the charge sector 
and the spin sector; depending on the magnitude of the hole- 
doping parameter, they exhibit both insulating and supercon- 
ducting behaviors in the charge sector, while in the spin sector 
it always displays short-range spin correlations^. 

The effects of mobile holes in the spin-gapped background 
are interesting in their own light not only in purely theoret- 
ical context^ but also in the experimental point of vie w^^'-^^ . 
However, the impact of mobile holes on the (hidden) topologi- 
cal properties has been little studied. In what follows, we will 
show that the sVBS states possess a kind of non-local topo- 
logical order in the spin sector as well as local superconduct- 



ing order in the charge sector, the latter of which is already 
known. While various (ordinar3/}correlation functions have 
been investigated akeady in ReflsTi dynamical properties, as 
exemplified by magnetic- (triplon) and charge (spinon-hole 
pair, specifically) excitations, are yet to be understood and will 
also be addressed in this work. 

In the sVBS models, supersymmetry (SUSY), i.e. rota- 
tional symmetry of boson and fermion, is realized as the sym- 
metry of bosonic spins and fermionic holes. Such SUSY of 
the sVBS states is exact regardless of the magnitude of hole- 
doping parameter, and their parent Hamiltonians can be read- 
ily constructed based on such (super)symmetry. Thus, the 
sVBS models enable us to systematically study hole-doped 
antiferromagnets on such a firm mathematical background. To 
this end, we develop a supersymmetric version of the matrix 
product state (MPS) representation of the VBS-type states—. 
Since the sVBS states generally contain fermionic degrees of 
freedom, we generalize the MPS formalism to include both 
fermionic and bosonic operators. This supersymmetric MPS 
(sMPS) representation is useful not only in the sense of com- 
putational efficiency, but also from the topological-order point 
of view as the emergent edge degrees of freedom, which char- 
acterize the topological features, are automatically incorpo- 
rated in the MPS formalis m'^'^^i-^^ . It should also be men- 
tioned that the MPS formalism, which has been introduced^ 
originally as a special class of quantum ground states with 
short-range correlations, is now believed to be a natural frame- 
work to represent entangled quantum many-body states in 
ID^ML. In a similar sense, the sMPS formalism would be 
applicable not only to the sVBS states to be investigated in 
this paper but to a wider class of entangled many-body states 
that contain fermionic degrees of freedom. 

This paper is structured as follows. In section |II] we intro- 
duce type I and type II sVBS states and summarize some ba- 
sic features. In section |III1 by including fermionic degrees of 
freedom, we develop the sMPS formalism, and apply it to the 
calculations of physical quantities of the type I sVBS states. 
The generalized hidden order is proposed and the string or- 
der parameter is evaluated in section |IV] In section |V] we 
calculate the gapped excitation spectra of the magnetic- and 
the charge (i.e. hole-pair) excitations on sVBS chains within 
the single-mode approximation. In section |Vl] we proceed to 
the analysis of type II sVBS states and derive the hole-doping 
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behavior of various physical quantities (e.g. superconducting 
order parameter and string correlation). The stability of the 
hidden 'topological' order found in these states is discussed 
from the point of view of the entanglement structure in sec- 
tion Vll. Section lVlIIl is devoted to summary and discussions. 



II. BASIC PROPERTIES 

Before proceeding to the detail analysis, we quickly review 
the basic features of the sVBS states in this section. 



A. Type I SUSY VBS states 

In what follows, we analyze two types of sVBS states. The 
first is the sVBS states with UOSp(l |2) supersymmetry"^- pro- 
posed recently in Reflsil (see Appendix |A] for a very brief 
summary of supersymmetry), which we shall call type I: 



In particular, for M = 1, each site can take two spin values 



5 -1 i 

— ^, ^, 



(5) 



and the local Hilbert space is spanned by the five (4M+1, in 
general) basis states 

|1) = ^af Ivac), |0) = al6,l|vac), = ^^iVac), 



|t)=al//|vac), U) =6://|vac) 



t/ti 



(6) 



Mathematically, these constitute an M—1 SUSY multiplet, 
and hence we use the name 'Type F. In addition to the lo- 
cal physical degrees of freedom on each site, the following 
emergent degrees of freedom localized around the edges {edge 
states) will play an important role: 



It)) I vac), i;)) =6t|vac), 



= /t|vac). (7) 



isvBs-i)=n(«i^^-^H-'^///i; 



vac) 



(1) 



where {ij) signifies a pair of adjacent sites (i, j) and r stands 
for the hole doping parameter The operators a^, bi and fi re- 
spectively are a pair of the standard Schwinger bosons satisfy- 
ing [fli, aj] = fej] = 5ij and a (spinless) fermion satisfying 



As we will see in section IIII Al the ground state of a finite 
open chain is 9-fold degenerate (corresponding to the 3x3 
matrix for the M = 1 type-I sVBS states). 

The M ~ 1 type-I sVBS chain interpolates between the 
two VBS states in the two extremal limits of the hole doping: 
at r — 7^ 0, |sVBS-I) reproduces the original spin-1 VBS state 
|VBS)i2 



{/j, /J} = 5^^. The vacuum |vac) is annihilated by both the |sVBS-I) |VBS) = JI^^-I^j+i ^ ^l«i+i)|vac), (8) 



boson and the fermion: a|vac) = 6|vac) = ./|vac) = 0. Since 
the fermions always appear in pairs of the form .//./J (i, j are 
adjacent), the sVBS states can be regarded as the hole-pair 
doped VBS states. One can easily see that the state |sVBS-I) 
is UOSp(l |2)-invariant from the invariance of the matrix 



while, in the limit r ^ oo, |sVBS-I) reduces to the Majumdar- 
Ghosh (MG) dimer stately |MG) 



7^, = 




|sVBS-I)^|n/t||MG), 



(9) 



(2) 



where |MG) is either of the two dimerized states of the MG 
model^: 



used to construct |sVBS-I) (The parameter r is absorbed in 
the renormalization of /. To see how the matrix TZi is related 
to |sVBS-I), see section IIII Al ). and hence |sVBS-I) has the 
UOSp(l|2) symmetry [See Appendix IA ll for more details]. 

The type-I sVBS states^' ([T]), that contain (fermionic) hole 
degrees of freedom as well as the (bosonic) spin ones, are a 
generalization of the standard spin-S* VBS statesii^ii . In the 
type-I SVB states ([l]), the total particle number at each site is 
conserved: 



|MG) 



vac) 



(10) 



For larger M, |MG) should be replaced with the inhomoge- 
neous VBS statesii where the number of valence bonds alter- 
nates from bond to bond. 

According to the spin-hole coherent state formalism^^, the 
sVBS state is expressed as 



zM = a\a, + b\b, + flu , 



(3) 



(ij) 



(11) 



where the lattice coordination number z is 2d for the d- 
dimensional hypercubic lattice (in what follows, z ~ 2 un- 
less otherwise stated). The integer zM plays a role of the spin 
quantum number 25 in the usual VBS states. Since /t/ takes 
either or 1, the following two eigenvalues are possible for 
the local spin quantum number Si = h{a\ai + b\bi): 



which is simply obtained by replacing the operators a,b, f 
with their corresponding classical counterparts u, v, 9. {u, v 
are Grassmann even quantities, while 9 is Grassmann odd.) 
From the Grassmann odd properties of 9, ^'svbs i can be 
rewritten as 



M, M 



(4) 



*sVBS4 = cxpl -Mr^ 



ViU, 



■ $ 



VBS 



(12) 



3 




+ 




FIG. 1. The type I sVBS is a superposed state of hole-pair doped 
VBS states. With finite hole-doping parameter r, all of the hole-pair 
doped VBS states are superposed to form the sVBS state, and the 
sVBS state exhibits the superconducting property. At r = 0, the 
sVBS state is reduced to the original VBS state (depicted as the first 
chain), while r oo, the sVBS state is reduced to the MG dimer 
state (depicted as the last two chains). 



where 'i'vBS = Y[{ij)i'^i''^j ~ ViUj)^ is the spin coherent 
state representation of the original VBS state. This expres- 
sion reminds the BCS wavefunction of the superconductivity; 

iBcs) = Uka+9kcici,)\o) - cMEkdkcicUm with 

electron operator Ck and coherence factor (See chapter 2-4 
in Ref. 1471). In both ^'svbs i and |BCS), the fermions alway 
appear in pairs and the wavefunctions can be expressed by a 
superposition of such fermion pairs, as demonstrated by ex- 
panding the exponential (See FiglT). 



B. Type II sVBS states 

The type II sVBS state is an extension of the previous series 
of VBS states (type I) and now contains doped (antisymmet- 
ric) bound pairs of two species of holes. The inclusion of 
two species of holes / and g allows us to write down a wave- 
function more symmetric with respect to the bosonic- and the 
fermionic degrees of freedom. Now, we introduce the type II 
sVBS states of the form: 

\sWBS-ll) ^l[{alb]-bla]-rfjg]-rgjfl)''\vac), (13) 
which is associated with another matrix: 



Tin = 






1 




































-1 






(14) 



the type II sVBS states have the UOSp(2|2) symmetry larger 
than UOSp(l |2) symmetry of the type I sVBS states. 

We have two species of fermions, and the total particle 
number at each site i is constrained by 



zM = a\a, + b\b, + f}f,+g\g. 



(15) 



where z is the lattice coordination number. Since the eigen- 
values of nf{i)=flfi and ng{i)=g\gi can take either or 
1, in the type II sVBS chain {z = 2), the following four 
eigenvalues are allowed for the local spin quantum number 

= Ua\a, + b%): 



S,=M, M-i A/-i 



(16) 



which respectively correspond to the possible combinations of 
the fermion numbers: 



infii),n,{i)) = (0,0), (1,0), (0,1), (1,1) 



(17) 



In particular, for the M ^ I sVBS chain (i.e. z = 2), the 
possible values read 



S^ 



1 1 

^' 2' 2'0 



(18) 



Therefore, the local Hilbert space is spanned by the following 
nine basis states 

|1) = 7|«lVac), |0) = a|fol|vac), j-l) = ^^flvac), 

|t)=al/;|vac), \i)=blfl\ve^c), 

|t') = ay|vac), |;')=foy|vac), 
tfti 



|0')=5|//|vac) 



(19) 



The name 'type IF is indicative of an JV~2 SUSY multiplet 
formed by these states. Again, |vac) denotes the vacuum with 
respect to (a, b, /, g). The edge states are now given by 



It)) = at I vac), i;)) =6t|vac), 
|0)) = /t|vac), |0'))=.gt|vac). 



(20) 



and, correspondingly, there appear 4x4 = 16 degenerate 
ground states for the i\/ = 1 type-II sVBS chain (see section 
[Vlfor the detail). 

The M = 1 sVBS chain has the following properties. As 
in the type I sVBS, it reproduces the pure spin VBS state for 
r -> 0: 



|sVBS-II) ^ |VBS) = Uialbl, 



^^al+i)|vac) 



(21) 



The new fermion 
relations {gt,gj} 
type II sVBS state 
due rescaling r 
the (spinless) hole, 
local sites fjgl\Q) 
type I sVBS states. 



gi satisfies the standard anti-commutation 

~ {fi^9j} — 0' 6tc. Apparently, the 
reduces to the type I after gi fi (and the 
ir). With inclusion of another species of 
. in the type II VBS states, there appear the 
with spin-0, which are not realized in the 
As we will show in the end of this section. 



On the other hand, when r — oo, it reduces to the totally 
uncorrelated fermionic (F) state filled with holes: 



|sVBS-II) ^ IF-VBS) EE n(//5j+i +5j/i^i)|vac) 

i 

= ±n//5l|vac) 



(22) 
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(the sign factor depends on both the parity of the system size 
and the edge states). Here we have assumed the open bound- 
ary condition^ and the normahzation of the edge states given 
in section IVI Al 

This properties are quite similar to those of the BCS state; 
at gk — > 0, the BCS state reduces to the electron vacuum (no 
fermions), while for gk — > oo, it coincides with the Fermi 
sphere (filled with electrons). In this sense, the type II sVBS 
states look more similar to the BCS state than the type I VBS 
states. As in the previous case, one can pursue the analogy to 
the BCS wave function by using the spin-hole coherent state 
representation of |sVBS-II): 



*sVBS-II ==11 C^i^J 
<i]> 



n( 

expf -Mr ^ 



ViUj — rUiTij 



- rrii0j) 



M 



<tj> 



■ expf -Mr^ ^ 



<tj> 



{u,Vj 



■$VBS. (23) 



Expanding the exponentials, one can easily see that with finite 
r, the type II sVBS states can be expressed as a superposition 
of the hole-pair-doped VBS states and that the system exhibits 
the superconducting property. However, unlike type I, type II 
sVBS states have no spin degrees of freedom at r — > oo. The 
intuitive picture of the M = 1 type II sVBS chain is depicted 
in Fig|2] 

Before concluding this section, we give a remark about the 
symmetry of the type II sVBS state. In Ref.lU an apparently 
different form of the sVBS states 



IsVBS') 



fy]-rglg]r\vac), (24) 



has been introduced. The state |sVBS') is manifestly invariant 
under the UOSp(2|2) transformation, since it is constructed by 
using the UOSp(2|2)-invariant matrix 





FIG. 2. (color online) The type II sVBS states are expressed as super- 
position of the hole-pair doped VBS states. Unlike the type 1 sVBS 
states, the spinless sites, depicted by the large white circles with dou- 
ble holes, appear. The MG states are realized in the "middle" of the 
sequence. The original VBS state and the hole- VBS state are respec- 
tively realized in the first and last lines. 



As the phase of the operators can be chosen arbitrarily, flip- 
ping the sign of them does not affect physics. Therefore, both 
|sVBS-II) and |sVBS') have the same symmetry UOSp(2|2) 
in common and are physically equivalent; all physical quanti- 
ties take completely identical values for these two states^. 



7^' 




(25) 



(r can be absorbed in the normalization of / and g.) In fact, 
the two sVBS states |sVBS-II) and |sVBS') are physically 
equivalent. By the unitary transformation 



V2U 1 



the fermion-pair part of |sVBS-II) is transformed to 

flgj + glfj^f^fj-glg]- 



(26) 



(27) 



Then, we flip the sign of either gj or gj to recover the correct 
form of the fermion-pair part in |sVBS'): 



flfl 



gig} 



flf] 



(28) 



III. SUSY- VBS STATE-I 

In the following sections, we consider the sVBS states de- 
fined on one-dimensional (ID) chain. A simplest SUSY- 
extension of the ID spin-1 AKLT (VBS) statei*^ is defined 
as (Af = 1, z = 2 in eq.O): 



(29) 



|sVBS-I) • • )(a]_i6] - 6]_ia] - r flj]) 

(a]fe]+i-fe]a]+i-r/t/t^,)(...)|vac) . 

The (non-hermitian^) parent Hamiltonian for the SUSY 
(UOSp(l|2)) VBS model is given as^i: 

ni = l SVBS = {^3/2^3/2(C,., + l) + V2P2{C,,j + l)] , 

(30) 

where Cj.j+i and P;(Cj.j+i) respectively denote the 
UOSp(l|2) Casimir operator on a two-site cluster + 1) 
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(see eqs.( IA3IA7l l for the definition of Casimir operators) and 
the projection operator onto ^tot = I subspace (note that the to- 
tal superspin Ztot of two I = 1 superspins can take all integer- 
and half-integer values between and 2; see eq. (IA8l l). For 
the positivity of the Hamiltonian, we require V3/2, V2 > 0. 
Specifically, the local Hamiltonian /ij.j+i is given by the fol- 
lowing fourth-order polynomial of the Casimir Cj.j+i: 



h{C) = 



3/2 



6 

V 45 



70,' 



2V2 
63 



3V2 
70 

315 



431A 



3/2 



90 
2V, 



3/2 



45 



(31) 



A. Matrix-product representation 

First let us briefly recapitulate the basic properties of a 
generic (bosonic) matrix-product state of the following form 
(see, for instance, Refs. [l5| and Hi] for recent reviews of the 
matrix-product representations): 



IMPS) = 



(32a) 



where the matrix Aj consists of state vectors at the site-j and 
its size is determined solely by the size of the auxiliary Hilbert 
space and is independent of the number of sites^. The state 
|MPS) in general is not normalized and we reserve the nota- 
tion |MPS) (and |sVBS)) for the unnormalized states. Ground 
states which can be expressed in this form may be generically 
expected to have finite degeneracy. For example, the ground 
state of the AKLT model, which is expressed by the spin-S" 
VBS state, is showi>i« to have (S'+l) x (S'+l)-fold degener- 
ate, when the model is defined on a finite open chain. When 
the system is defined on a periodic chain, we have to take the 
trace over the matrix indices: 



IMPS) 



PBC 



Tr 



(32b) 



Below, we shall see that the expression eq.( l32bl i should be 
modified when A contains both bosonic degrees of freedom 
and fermionic ones. 

Now let us construct the matrix-product representatio n^^'^^ 
of the type I (UOSp(l|2)) VBS state When 
the Schwinger-boson/fermion representation of the state is 
known, the simplest way'' would be to find an operator- 
valued matrix in such a way that everytime when we multiply 
a new matrix (say, gj+i) from the right the (SUSY) valence- 
bond operator 



is inserted between the previous right edge (site-j) and the 
newly added site (j+1). To this end, let us introduce the 
'spinor': 



t^t 



in terms of which the above UOSp(l |2) valence bond can be 
written compactly as: 



(34) 



('t' denotes the transposition). The 'metric' TZi has been de- 
fined as 




Then the sVBS state 

(a,/3 = 1,2,3): 



7^I = I -1 I . (35) 
is written as a string of 3 x 3 matrices 



L-l 



|sVBS)„^ = {nii;ir Y[mniij,+i)ijl\yac) 



= {Tl,^iri^{ ■ n ^"/^'V.': • (7^IVLV'^)|vac) 



(^1^2••■^L)a^ 



(36) 



where 

A. : 



7^I-0J • V'j|vac)j 







|o), V21-1), 

= I -^/2|l), -|0), -^It), 



(37) 







= ^ T^^\m)\m)+ r(P)(TO)|m) . 

m= — 1,0,1 rn=t,i 

The 3x3 matrices F*^^^ and F'^^ respectively denote the 
bosonic- and the fermionic part. The edge operators TZiijji = 
(6| , —a\, —y/rfly and tpL appearing respectively on the left- 
and the right edge represent the three possible edge states 
(spin-up/down and hole) on each edge. 
Following the same steps as the above for 



(sVBS-I| ==(vac|(- • • ){ajbj+i - bjaj+i - r fj+Jj) 
{aj-ibj ~ bj-ittj - r fjfj-i){- ■■ ) , 



we obtain 



„^(sVBS-I| 



with 



A] 



(vac|7/'*7/']7^' 



t-nt 



ajbj 



(vac 



■ aIa\ 



Pa 



(38) 



(39) 




(40) 



(33) where ijj* = {aj,b.j,y/rfjy. 
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By construction, it is obvious that all the nine matrix ele- 
ments of the following string of A-matrices: 



blal 




(41) 



are the (zero-energy) ground states of the parent Hamiltonian 
^j-j+i- That is, the product (S'jLi gives the the 
ground states of the M = 1 sVBS model on an open chain 
with length L. Here it is important to note that we are free 
to choose the polynomials ((6| —a\ —y^/j) from the left 
edge and {^a\ b\ y/rfj^) from the right) appearing at the 
edges. As will be discussed in section ITlIB I this leads to a 
remarkable feature of the VBS-like systems-ec/ge states. 

In constructing the sVBS state on a periodic chain, one 
has to treat the fermion sign carefully and one sees that the 
trace operation used in the standard MPS representation (I32bb 
should be replaced with the supertrace (see AppendixlBl): 




IsVBS) 

periodic " 

where the supertrace here is defined as 

STt{M) = Mil + M22 - M33 



(42a) 



(42b) 



From these A-matrices, we can calculate the following 9x9 
T-matrices (transfer matrix): 

T{a,a-J,/3)=A*{aJ)A{a,f3) 





/I 











2 








r\ 







-1 








































— r 
















-1 
















2 











1 








r 



















r 






















-r 



















r 



















\r 











r 










(a, a 






1,2,3) , 









(43) 



where A* is obtained from Ahy \-] 
jugation. The eigenvalues of T are 



-I(x3),-ir(x2),ir(x2), 



and complex con- 



1 



V8r2 + 9 



(44) 



any finite r, by 



i (3 + VSr^ + 9) 



(45) 



In the limit r — > 00, another eigenvalue (3 — VSr^ + 9) /2 
becomes degenerate with the above. 

The use of the supertrace in eg. ( 142 at modifies the expres- 
sion ( |B8t of the norm for the periodic system to: 



(46a) 



(MPS|MPS)pBc = 5]sgn(a)sgn(/3) 



where 



sgn(a) = 



1 for a = 1,2 
— 1 for a ~ 3 . 



(46b) 




FIG. 3. (Color online) Plot of absolute values of the five different 
eigenvalues of T. The largest eigenvalue is always unique and non- 
degenerate. 



B. Edge states 

Now we would like to mention an important feature of 
the VBS-like states defined on an open chain. From the 
expression i4T[ . it is clear that the nine degenerate ground 
states correspond to different choices of the edge polynomi- 



als {bl, 



-Vrfi) and {al, bl, ^/rfl). In fact, we can 



explicitly indicate the edge-dependence of the ground states 
as follows: 



|sVBS)open 




(47) 



The largest eigenvalue which is relevant in determining the 
physical quantities in the thermodynamic limit is given, for 



From this, we can readily see that the matrix indices of the 
MPS are directly related to the edge states. It is instructive to 
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calculate for various edge states |sVBS)ipen^''\ InFig.H 
we plot the local magnetization for three left edge states 
Sl (with the right edge state Sr fixed). 



MPS formulation, this is a direct consequence of the fact 

In fact, this property greatly simplifies the calculations below. 
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FIG. 4. (Color online) Plot of (S-) (r = 0.3) for various (left) edge 
states 'hole', 'f and 'J,' (with the right edge state fixed to sr —f). 
The system is non-magnetic in the bulk and magnetic moment exists 
only around the edges of the chain. 



C. Spin-spin correlation 

Now that we have obtained all the necessary matrices, we 
can follow the steps described in section lB^ to calculate vari- 
ous correlation functions. 

The ordinary spin-spin correlation function (>S'"S'"_|_„) 
reads: 



2 (r2 + 3 + \/8r2 + 9) 



V8r2 + 9 (3 + V8r2 ^ gj 
13r2 + 24 + (r2 + g) \/8r 



(for n = 0) 



(49a) 



2V8r2 + 9 (3 + V8r2 _^ g) y 3 + V&r^ + 9 
(for n > 0) . 



(49b) 



The exponentially decaying factor defines the correlation 
length^: 



A remark is in order here. One may think of the above edge 
moments (s = 1/2 moment or a hole) as independent physical 
objects and conclude that the (SUSY) VBS states are orthog- 
onal with respect to these edge states. However, this is not 
true; in fact, the above edge moments are emergent objects and 
sVBS states with different edge states have finite overlaps with 
each other, which are exponentially decreasing as the system 
size L. That is, two VBS states with different edge states are 
orthogonal to each other only in the infinite-size limit. In the 

I 



CspinC?-) ^ = log ■ 



9 



(50) 



which is monotonically decreasing in r. In the pure AKLT- 
limit r ^ 0, it reduces to the well-known resultsi*^: 



- (—Y 

3 V 3 / 



for n ~ 
for n > 



For later convenience, we calculate the static structure factor S°'"{k). The result is given as: 



4 (2r'' + 17r2 + {3r^ + 6) V&r^ + 9 + 18) (1 - cos(fc)) 



V8PT9 (V8r2 ^ 9 + 3) |4^2 3^8r2 + 9 +11 + 2 (V8r2 + g _^ 3) cos(fc)} ' 



(51) 



(52) 



D. superconducting correlation 



ber or 1 at the site fc; see Fig|5]i: 



In order to handle the operators containing fermions, we 
have to generalize the general recipe presented in Appendix 
IbI Take for example the hole-pair creation operator—: 



(53) 



i-lf'A, . (54) 



(-1)^^-M,_i and/. 



Next a string (-1)^^ Ai (g) • 
interchanged and this multiplies the matrices Ai, Aj^i 
additional (—l)^*" -factors to remove the fermion sign except 
at the site j. Therefore, we need four more matrices 



In order to apply the method presented in sections lBT] and lB2l T'^f'' {^a a- jS /3)^A*{a i3){Op)A{a (3) 
„ f A ™„t,; — A ^. -^A - has to be moved to the ' ' ' ' 1 ' 



first a string of A-matrices A 



left of /J^j^ and through this procedure it acquires a Jordan- 
Wigner-like phase 01=1 (^1)^'" (-^k counts the fermion num- 



T<^f' (a, a; ^, /S) ^ A* (a, p) {Op {-if} A{a, (3) (55) 
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By using these, the numerator of (Aj) is calculated as: 
Also interesting are the hole density 

("hole) = (/;/,) 

and the hole-number fluctuation 



6nm, = ^{f]h)-{f]hY 



(56) 



(57) 



(58) 



By using the method described above, we can readily calcu- 
late these quantities. For instance, the hole density in the bulk 
system is computed as: 



r2(5 + \/8r2 + 9) 



8r2 + 9 + (r2 + 3) VSr^ + 9 



(59) 



As is clearly seen in the inset of Fig. |6] near the edges of an 
open chain, the hole density is different from the bulk value 
and approaches exponentially with the 'healing length' given 
by 



iho\e{r) ^ = log 



9 + 3 



(60) 



Note that this is different from the spin correlation length 
Cspin(?') in eq.dSOl) and the superconducting correlation length 



6c(r) -1/log 



9 + 3 



2r 



(61) 



defined by the exponential decay of the singlet off-diagonal 
ior 

/] IMPS) 



correlation functior>^ {{a^bj+n ~ bjaj+n)f]f]+n)- 



(-If (-lY 



j-i 

(a) 



j+i 



(b) 



FIG. 5. Action of fermion operator on the MPS. (a): Due to the 
fermionic anticommutation relation, extra factors ( — 1)^ appear in 
the ^-matrices on the left of site-j. Accordingly, a new transfer ma- 
trix (c) is necessary as well as the standard one (b) when we calculate 
expectation values containing fermionic operators. 



In Fig|6l we plot the expectation value of the hole-pair cre- 
ation operator: 



(62) 



together with the hole density (rihoie) 



and the 



From 



hole-number fluctuation (5nhoie = \/ ifjfj) ^ if] fj) 
r = (5 = 1 VBS limit) to r ^ oo (S" = 1/2 Majumdar- 
Ghosh limit), the hole density is monotonically increasing. 
When r = and r — > oo, the hole number fluctuation is 
suppressed (rihoie takes definite values and 1, respectively) 
and consequently the 'superconducting correlation' becomes 
zero. This is consistent with what we expect from the analogy 
to the BCS wave function pointed out in Reflsil 




0.4 „ 0.6 



1 + r2 



FIG. 6. (Color online) Plot of 0,c = (A^), the hole density 
(?ihoic(j)) = {.fjfi) 'irid the hole-number fluctuation (fjfj) — 
ifjfj}^ as a function of r. Here the bulk values are plotted. In- 
set: Profile of the hole density (r — 0.5) for a finite system (L — 20) 
with different left edge states (t, and 'hole'). Only the left edge 
state is changed with the right one fixed to sr The hole density 
approaches exponentially to the bulk value as we move away from 
the edge. 



IV. HIDDEN ORDER 

A. Generalized Hidden Order in sVBS states 

The hidden order is a generalized concept of the Neel order 
For 5 = 1 antiferromagnetic spin chain, the Neel order is 
depicted as 

••• + - + - + - + - + ••• (63) 

Here, + stands for Sz = +1, and — for = —1. In the 
sequence, + and — are alternating, representing the classi- 
cal antiferromagnets. A typical Sz sequence of VBS chain is 
given by 







00 







(64) 



When we remove zeros in the sequence, we arrive at the 
usual Neel order This is the hidden (string) order observed 
in gapped antiferromagnetic spin liquids-lsS^. The hidden or- 
der is a non-local order, since the removing zeros is a global 
procedure. Since in the sVBS states one-hole states carry one- 
half spins at each site, Sz — 1/2 and —1/2 generally appear in 
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the sequence. The locations of such one-half-spins are, how- 
ever, not completely random; The following procedure reveal 
the existence of a generalized hidden order in the s VBS states. 
A typical Sz sequence of sVBS states is given by 

■ ' ■ ^ tj' ^ ^ + -00 + ij, ■■■ (65) 

First, we search the spin-half sites from the left and whenever 
we encounter a pair of spin-half sites we sum the two S,- 
values to replace the pair with a single site having the effective 

5,(e.g. ;;^-): 



+ 00 - + 



+ 







(66) 



Then, we remove the zeros in the sequence to obtain the stan- 
dard Neel pattern: 

••• + - + - + -••• (67) 

This argument leads us to conclude the existence of (gen- 
eralized) hidden order in the sVBS states. By the SU(2)- 
invariance of the sVBS state, the same is true for the Sx- 
sequence as well. The hidden order is "measured" by the non- 
local string order parameter^'*. In sections H V CI and [VI D| we 
explicitly calculate the string order for the type I and the type 
II sVBS states, respectively. 



B. Matrix-product representation and hidden order 

Before proceeding to the actual calculation of the string 
correlation, we delve the hidden order inherent in the sVBS 
state from the MPS point of view. Since the condition for the 
string correlators to have finite values is known in a general 
and mathematical manner—, we give here a more physical ar- 
gument. 

To clarify this hidden structure in the spin configuration, 
let us pick up an arbitrary site j and consider the partial sum 
of S^s contained in the block between the left edge and the 
site-j: 



To see what (fTOl i implies, it is suggesting to plot S^^^lj) as a 
sequence of steps. Namely, we assign a local height variable 
hj = 5,ot(j) to a bond to the right of the site j. Then, the 
local spin value Sj is expressed as a step hj — /ij-i between 
the adjacent heights. It is obvious that this height plot is in 
one-to-one correspondence to the original {S^} configuration. 
Eq.dTOli shows a set of possible heights (i.e. 5tot(j )) ^ given 
site j. For instance, if the left edge state is t, the correspond- 
ing states are contained in the first row of (iTOt and one readily 
sees that only 0, 1 and 1/2 are allowed for the sVBS state. 
Figl?] shows a typical height configuration corresponding to 
the usual VBS stated (a) and its SUSY counterpart (b). Strik- 
ingly, the height configuration is always meandering between 
the height-0 and the height- 1 (although the absolute height of 
the meandering line depends on the left edge states, the height 
configuration is always confined within a region of width 1). 
The same reasoning applies to the general spin-S* VBS cases 
and we can show^i^ that the height configurations are confined 
within a region of width S. This is highly non-trivial since in 
the ferromagnetic state we have an ever going-up steps. This 
'almost flat' feature of the VBS state has been first realized by 
den Nijs and Rommelse^"* for the 5 = 1 case. 

In the case of 5 = 1, one can strengthen this statement; in 
any spin (or height) configurations satisfying the above prop- 
erty, S*^ = 1 and —1 occur in an alternating manner when 
the intervening Os are neglected (see Fig|7ja)). This may be 
viewed as a diluted Neel order In the standard Neel state, 
we can insert an alternating phase (—1)^"' to make the cor- 
relation between the two spins and Sj ferromagnetic. In 
the diluted case, on the other hand, we can easily see that the 
string operator Y[k=i exp(i7r5^) will do the job and that one 
can use the following order parameter (string order parame- 
ter) detects the Haldane stata^i^: 



/ J+ri-l \ 

0-„g = jim(^^| n cxp{»7r5,n5|+„ 



(71) 



k=l 



(68) 



In considering the possible values of it is convenient 

to consider the MPS for the block: 



{^1 



(69) 



Since the sVBS state on any finite subsystem ( 1361 ) is made up 
of a product of (SUSY) valence bonds ( |34] | carrying 5" = 0, 
the above S^^^{j) is determined only by the edge states of the 
subsystem 



{^1 



15= 



,(J)=0) 



\S^oAj)=-1) |5f,(j>-l/2)^ 
|5.^ot(j)=0) \S^,ij)=l/2) 



|5.^o,(j)=l> 

|5^-„.(j>l/2) \SUj)= - 1/2) \S^Jj)=0) 



(70) 



For the spin-1 VBS state, it is evaluated^ exactly as (2/3)^ 
('2/3' comes from the probability of having non-zero S^). 

In the SUSY case, the situation is slightly more compli- 
cated since we have height- l/2s corresponding to sites with 
one hole. However, if we note that the holes appear always in 
pairs, we can easily see that the insertion of hole-pairs (which 
carry = 1/2) does not affect the string part 



j+n-l 



j+n-1 



Yl cxp [ittSI] = exp lii: S'k 



k=3 



and we may expect that string order persists in the SUSY case 
(r 7^ 0) as well (see Fig|7jb)). 



10 



(a) 



' k=l 


- 


1 








- 


1 


- 


1 
















Sl = 





(b) +, 

+1/2 




1 -1/2 -1/2 1 -1/2 1/2 -1 ! -1 



hole pair hole pair 



FIG. 7. (Color online) Height plot of typical spin configurations in 
spin- 1 VBS state (a) and M = 1 sVBS state (b). Note that heights 
are confined within a region of width 1. Although simple 'diluted' 
Neel picture does not hold because of the presence of hole pairs, still 
we can find string order when hole pairs are grouped together in (b). 



C. String correlation 

The finite-distance string correlation function^ 



CstringO ; n) 



Sj exp 



j+n-l 

- E 

k=j 



s. 



j+n 



(72) 



observation that the spin-1 Haldane state is adiabatically con- 
nected to the spin- 1/2 dimer stated. Meanwhile, the type II 
sVBS chain {M = 1) is reduced to the hole- VBS chain with 
no spin degree of freedom at r — > oo, and hence the string 
order vanishes completely in this limit. 



In Ref.m a SUSY-analogue of the higher-S* VBS states is 
discussed as well. The ordinary spin-5 VBS states obtained in 
the zero hole-density (r 0) limit are known to exhibit dif- 
ferent topological properties according to the parity of spin- 
S; the string order parameter vanishes for the even-spin VBS 
states while it is finite for odd-SiSil. In this sense, it would 
be interesting to calculate the string order parameter for 
the generalized sVBS states. As is seen in eq.([T]i, the role of 
spin 5* is played by an integer M (superspin) in the SUSY 
case. For all M, we can construct the matrix-product repre- 
sentation of the AZ-sVBS state by using {2M + 1) x {2M + 1) 
matrices (see Appendix |C]l and after straightforward evalua- 
tion we obtain the results shown in Fig. |9] As is expected 
from the previous studies, the r = value of 0^^^^„ vanishes 
for even- A/. When the hole pairs are doped, on the other hand, 
the string order revives. In section [Viri we will interpret this 
from the point of view of symmetry-protected topological or- 
der 



can be evaluated in a similar manner In the case of open 
chains, it suffers from the boundary effects. However, if we 
consider the case where both end points j and j + n aie in- 
finitely far from the chain edges, the expression simplifies a 
lot. In general, it contains exponentially decaying parts 



V8^ 



9-3 



9 + 3 



(73) 



as well as the constant (i.e. 
Fig.©: 



long-range-ordered) one (see 



stnng 



(r) 



4 {r" -f- 14r2 + 18 -f 2 (r^ 4- 3) 

(8r2 -h9) (V8r2 -f 9 + 3)^ 



(74) 

Only in the limit r ^ 0, the exponentially decaying parts 
disappear and the string coiTelation function becomes con- 
stant 4/9 (perfect string correlation). Note that the correla- 
tion length String is different from that (4pin) for the spin-spin 
correlation. With increase of the hole-doping parameter r, 
the effective spin magnitude gets reduced by the increase of 
the spin-half sites and accordingly the string order parameter 
monotonically decreases (see Fig.|8]l. 

At r — > oo, the type I sVBS chain {AI = 1) realizes the 
Majumdar-Ghosh dimer states with one-half spin degrees of 
freedom at each site and the string order parameter 0"^^^ 
reaches its finite minimum 1/16, which implies that the string 
order survives even in the r oo limit. This agrees with the 



e 


2.5 
f 2.0 




J 1.5 

5 


to 

■S 02 


1- 0.5 











r=0 





FIG. 8. (Color online) String correlation function in the bulk 
Cstring(cxD; n) for various values of r: (i) r = (top; pure spin 
AKLT), (ii) r = 5.0 (middle) and (iii) r = 20.0 (bottom). Note 
that for the pure spin AKLT model (r — 0), the string correlation 
function is constant 4/9. For r 7^ 0, the string correlation functions 
exponentially approach to the limiting values shown by dashed lines. 
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1 + r2 

FIG. 9. (Color online) The infinite-distance limit of the string corre- 
lation function O^rinB ~ liinn->oo Cstring(oo; n) for several values of 
M plotted as a function of r. Note that C'Sing('"=0) = for even-M 
corresponding to the vanishing of string order parameter for even-S. 



V. SINGLE-MODE APPROXIMATION 

In this section, we consider the dynamical quantities, i.e. 
low-lying excitation spectra by using single-mode approxi- 
mation. As is easily verified, the so-called Lieb-Schultz- 
Mattis twist^, which provides a basic picture of gapless low- 
lying excitations in half-odd-integer spin chains, does not 
work in the usual VBS state- . Instead, an excited triplet 
bond {crackion- a 'crack' in a solid of valence bonds) in 
the valence-bond solid gives, to good approximation, a phys- 
ical low-lying excitation. As has been shown by Fath and 
Solyom^'-, the crackions are equivalent to the triplon excita- 
tions created by spin operators S(k). 



A. Spin excitations 

Let us start by investigating the action of local spin opera- 
tors 



Sj — cij bj , Sj 



b]a,, S]^l{a]a,^b]b,) (75) 



on the sVBS state. A little algebra shows that these spin oper- 
ators create triplet bonds around the site j (see FigfTOll: 



5+|sVBS-I) = |7^jL\>-|Vf^) 



S-flsVES-I) = 



>4{-i^j^>+i^r>} 



(76a) 
(76b) 



where IV'^^') and are obtained by replacing the SUSY 

valence bond (a]a]_|_j — — r by triplet bonds 



i^a^^i and (a|5j_(_^ +5|a]_|_^), respectively. This imphes that 
the triplon-crackion equivalence holds in the sVBS case as 
well. 



The single-mode approximation to the magnetic excitations 
is given by 



"SMA 



(fc) = - 



1 (sVBS-I| [[n,S°'{k)],S°'{-k)] |sVBS-I) 



2 (sVBS-I|S'"(fc)S'"(-fc)|sVBS-I) 
_ (sVBS-I|5'"(fc)-H g"(-fc)|sVBS-I) 
^ (sVBS-I|S'"(fc)S'"(-fc)|sVBS-I) 



(77) 



By the SU(2) symmetry, it suffices to evaluate wsma only 
for a = z and the spin index a will be suppressed here- 
after Using eq.( l76bb . the denominator (static structure factor) 
(sVBS-I|5"(fc)S'"(-fc)|sVBS-I) is evaluated as: 



(sVBS-I|^"(fc)5"(-fc)|sVBS-I) (= S'^k)) 

= i(l -cos/fc)(?A(°)('t)lV'^°H^)) 
where (k)) denotes the Fourier transform 



(78) 



\i,(^)ik)) = ^Yl 



-^fcr|^(0) 



Similarly, the local property of the sVBS states 
(sVBS-I|/i,j+i = /i,j+i|sVBS-I) = (Vj) 
implies that only the diagonal part survives: 

From this, one deduces: 

(sVBS-I|S'"(fc)-H S'"(-/c)|sVBS-I) 



(79) 



(80) 



1 



= -(l-cosfc)(V'f 
= i(l-cosfc)(V'(°)(fc)|H|^'W(fc)) . 
Eqs.dTSl) and dSTl i are combined to give 

„ (^(")(fc)|H|^(°)(fc)) 



/,(0)\ 



(81) 



(V'(o)(fc)IV'(°H^)) 

1(1 — cos k) 



(82) 



2 S^'{k) 



At this point, one may note a peculiar feature of the VBS- 
like states. Normally, a local excitation created by physical 
operators (e.g. Sf) propagates on a lattice by using the off- 
diagonal matrix elements: 



(83) 



In the VBS-like models, on the other hand, (^'1°'' l^lV-'j'^'') is 
diagonal by construction (all the diagonal elements are given 



/,(o)\ 
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by (^j^-* l^j'^'')) and excitations cannot use this chan- 
nel. Rather the non-trivial fc-dependence of ujsMA{k) comes 
only from the non-trivial overlap between the crackion states: 
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{(V'"'(fc)|V'"'(fc))} 



(1 — COS k) 



(84a) 
(84b) 



An important conclusion can be drawn from eq.(l82t: the phys- 

(fc)) becomes 



^SMA 



ical triplon excitation energy w^(fc)(< 
zero (i.e. gapless) as fc — > unless the static structure factor 
S''''{k) behaves like (fc 0). For any spin-S" VBS states 
and the sVBS states, we have checked that S^^{k) contains a 
factor (1 — cos fc) ~ fc^, which opens a gap at fc = 0. 



SUSY valence bond 



aici ri_c uui \u - 



(a) 



3^1 

iripleL bijr 



SUSY valence bond 



3[cnLt: uunu _ 



3-1 

_ spinon-hole pair i 

=f {- -e^^k^ + -B--€x>" } 

FIG. 10. (Color online) Action of local spin operator (a) and 
fermionic generator Ki (b) onto the sVBS state. The local operators 
Sj (a = x,y,z) and 1^^1,2 respectively create a triplet bond and 
a spinon-hole pair (crackion) on either of the two adjacent bonds 
U - and + 1). 



bond + 1) with a spinon-hole pair (aj/| + fjoj^i) (see 
Fig. [TOl b)). The excited state A'2|sVBS) is defined similarly 
with in the above expression replaced with . Then, the 
SMA excitation energy is given by an expression similar to 
eq.i 



Wsma(^) 



(sVBS-l|ifi(-fc)t-Hifi(-fc)|sVBS-I) 
(sVBS-l|Xi(-fc)tXi(-fc)|sVBS-I) 
(^(V2)0-)|H|V,(V^)(j)) 

(1/^(1/2) (fc) 1^(1/2) (fc)) ■ 



C. Fixing parent Hamiltonian 



(87) 



Before calculating the SMA spectra (|82] | and dSTl i. we have 
to fix the form of the parent Hamiltonian. As has been men- 
tioned in sectionHni the non-hermitian parent Hamiltonian for 
the SUSY (UOSp(l|2)) VBS model is given^by eq.i 



n, 



= lsVBS - ^ {^3/2-P3/2(Cj,j + l) + V2P2{Cj,j + i)} 



with the coupling constants I3/2, V2 > positive. 

The above form is not very convenient since it breaks her- 
miticity necessary for eq.(|79t and one still has one free pa- 
rameter even after the overall energy scale is fixed^' . Instead 
of using 'Hl=i sVbs, one may adopt 



n 



L=l sVBS 



n 



L = l sVBs'^i-=l sVBS 



(88) 



as the hermitian Hamiltonian^. One way to fix the remain- 
ing coupling is to require that the SUSY parent Hamiltonian 
should reduce in the r — > 00 to the standard (SU(2)) VBS 
Hamiltonian^ 



n 



5=1 VBS 



(89) 



B. Hole excitations 

A similar analysis can be done for the charged (hole) ex- 
citations which are always accompanied by spinon-like (i.e. 
S = 1/2) objects. These excitations are created by applying 
the two fermionic generators of UOSp(l |2) 



K2U) = lix-'f,b] - x/]a,) {x EE Vf) 



(85) 



to the VBS ground state. By using the explicit form of the 
ground-state wavefunction, it is easy to show 

Xi(j)lsVBS-I) = ^ {|vf_/i'^> - l^f Z'^)} , (86) 

where the crackion state is obtained by replacing the 

SUSY valence bond (a]fe]+i - b]a]^^ - r/j/j+i) on the 



However, this still has a problem; since some of the matrix 
elements in the fermionic sector have a factor 1 /r, the limit 
r — ?> 00 is divergent. Fortunately, this is not so serious. If 
we note that the ground states contain no fermion in the r ^■ 
00 limit, the most natural way is to require that the SUSY 
parent Hamiltonian projected onto the bosonic sector should 
coincide with the spin-1 VBS Hamiltonian ( f89t . This fixes the 
two coupling constants as^: 



3/2 



tanh r , V2 = V2 . 



(90) 



The spin-excitation ('crackion') spectrum obtained by using 
(IHl and dinil is shown in Fig HI] At r = (AKLT-limit), the 
dispersion reduces to the well-known ono^: 



'^sma(^') = ^(5 + 3cosfc) 



(91) 



For r 00, on the other hand, the spin excitation loses the 
dispersion. This is easily understood since the ground-state in 
this limit reduces to the translationally invariant combination 
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of two Majumdar-Ghosh states (see Fig. [U and the overlap 
between crackion states, which gives the dispersion of the spin 
excitations, trivializes (see ( |82] | and (I84ab ): 



(V'f ') « k] , = const. . (92) 



The charge excitation spectrum is calculated similarly by 
using eq.(l87]i. The result is shown in Fig.[T2] For r = 0, the 
spectrum is given by 



'^^^^''^ " 3(2- cos fc) • 




FIG. 12. (Color online) The excitation spectrum WsmaC^) of 
spinon-hole pair obtained by single-mode approximation (eq.l|87|)). 
This spinon-hole pair state is created by fermionic generator Ki ex- 
cept at r = 0, where the transition matrix elements of A'l from the 
ground state vanish. 



VI. SUSY-VBS STATE II 



A remark is in order here about the existence of the two 
different spectra aj^(fc) and uj^{k). One may naively expect 
a;*(fc) = uj^{k) as the supersymmetry relates the bosonic gen- 
erators S and the fermionic ones Ka- However, this relies 
on the existence of a 'unitary' transformation which linearly 
transforms the set of the SUSY generators onto themselves 
(adjoint representation). Since no such transformation exists 
here, we generally expect different spectra for the spin- and 
the charge sector as has been shown above. 



1.0 




FIG. 11. (Color online) The spin excitation (triplon) spectrum 
cJsma(^) obtained by single-mode approximation (SMA). At r = 
0, it reduces to the well-known dispersion ujsMAik) = 10(5 + 
3cosfc)/27 of the spin-1 VBS model. When r oo (Majumdar- 
Ghosh limit), dispersion becomes flat. 



Now let us add one more fermion species and consider yet 
another SUSY-VBS wavefunction which now includes two 
holes / and g. As has been mentioned in section HlBl the state 
contains two (spin) bosons (a, 6) and two fermions (/, g), and 
we may expect it to exhibit clearer spin-charge symmetry with 
respect to r = 1. 

The second generalized sVBS wavefunction (the case AI = 
1 of eq.(fT3])) is defined by: 

|sVBS-II) 

^{■■■){o]-A-^-A-^(fU9}+9Ufh} 

(94) 

As we have seen in section HlBl this state is based on the alge- 
bra UOSp(2|2) and one can construct the parent Hamiltonian 
in a similar manner to the type 1 case (based on UOSp(l |2)) 
(we do not give the explicit form here. The interested readers 
may refer the online supplementary material^.). 



A. Matrix-product representation 

We follow the same steps as in section |III] with a different 
metric matrix 



= 






1 




































-1 
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and the spinor (aj , b'- , ^/rfj , v^.g])' to obtain the MPS repre- 
sentation for the second sVBS state: 



/ 














[a] b] V-rf] V-rg]) 


|vac)j- 


V 










/ 


«]^' 




[b]f ^b]f] 










-a]b] ~^a]f] 
-^9]b] -rg]f] 








t 





V 




\ 

3 


-^rf]b] 


-^f]9l ) 



|vac)j 
(96a) 
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(bj ~a.j - 










( Oi^i 

(&.)^ 




-Vrajgj ^Vra^fj 


-ajbj 






-rfjgj 


Wrgjbj 


-VrgjCLj 


-rgjfj 



(96b) 

As in the first sVBS state, the supertrace is necessary for 
the periodic system: 



|sVBS-II) = STr 




(97) 



where STr(X) = A^n +7W22 -^33 TheT-matrix 
is a 16x16 matrix and has seven different eigenvalues Xi (see 
Fig.lll: 

{AJ ={-l(x3),-M'(x4),+ir(x4),-r2(x2),r^ 

1 1 (98) 



where /(r) = V r-^ + lOr^ + 9. Regardless of the value of r, 
the eigenvalue with largest modulus is: 



fir)) 



(99) 



Since the set of eigenvalues is invariant under r ^ —r, we 
can restrict ourselves to r > 0. 



B. spin-spin correlation 

Let us begin with the spin-spin correlation function. By 
using the method described in Appendix IB 21 it is straightfor- 
ward to calculate the correlation function (S'"S'"^„): 



fir] 



r) ( 2 \ 

\ r2+3+/(r) ) 



for n = 
for n > 



(100) 



FIG. 13. (color online) Plot of absolute values |Ai| of the seven dif- 
ferent eigenvalues of G. Since |Ai|s are symmetric with respect to 
r I— >■ — r, only the r > part is shown. 



In obtaining these expressions, it has been assumed that 
both end points (x and a; + n) are infinitely far from the edges 
(otherwise there will be another decaying factor coming from 
the edge effects). From these, we can read off the spin-spin 
correlation length: 



?spmM = l/log{(r2+3 + /(r)) 12] 



(101) 



which monotonically decreases from l/ln(3) (r = 0) to 
(r 00). 

The existence of the edge states may be best illustrated by 
plotting the local magnetization {Sj). 
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FIG. 14. (color online) Plot of local magnetization profile {S^) for 
different (left) edge states (with right edge fixed), (inset) Spin corre- 
lation length ^spin(T') as a function of r. It monotonically decreases 
as r is increased and approaches to zero like ^spin ~ 1/ log (l + r^) . 



C. Superconducting correlation 

Since the type-II sVBS state (O contains hole pairs on 
adjacent sites, we may expect that the pair amplitudes take 
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finite expectation values. As in section HIIDI we may define 
the following order parameters on general grounds: 



A// 


= {ajbj+i 


- feiaj+i)/]/J+i 




(102a) 


Af 


= {ajbj+i 






(102b) 


Af 








(102c) 



However, the first two are identically zero by construction of 
I sVBS-lI) . The only non- vanishing superconducting order pa- 
rameter 



(A 



fg\ 



(103) 



is plotted in Fig. [15] for various values of r. Also plotted are 
the hole (/ and g) number {nf,g) and the hole-number fluctu- 
ation (5nhole- 



(Jnhole = (?^hole> - {nho\e)^ ("-hole = rif + Ug) 



(104) 



The superconducting order parameter Osc is maximal at r w 
1.05 (or, rV(l + r2)~0.52). 
The superconducting correlation (hole-hole correlation) 

C//(n) ^ (a,&,+„ - 6,a,+„)(/]5]^„ + g]j]^J (105) 

decays exponentially with the correlation length 



Csc(f) = log 



r2 + 3 



2r 



(106) 




FIG. 16. (Color online) Plot of hole correlation Cl^{n) = 
(ajf)j+„-&3aj+„)(/Jg]_|_,j+g]/J_|_„) for various r. Due to thefom 
of the wave function, hole correlation identically vanishes when the 
distance n is even. Inset: coiTelation length ^sc('') of the hole corre- 
lation. 



D. String correlation 



Then, we proceed to the string correlation function. As 
in the previous case (type 1 sVBS), the string correlation ex- 
plicitly depends on the distance between the two end points 
through the exponentially decaying factor: 




1 +r2 



FIG. 15. (Color online) Plot of Osc = (A^), the hole density 
(nhokij)) = ifjfj) and the hole-number fluctuation (/J/j) — 
ifjfj)^ as a function of r. Here the bulk values are plotted. (In- 
set): Profile of the hole density (r = 0.5 or r^/(l + r^) = 0.2) for a 
finite system (L = 20). Only the left edge state is changed with the 
right one fixed to sr The hole density approaches exponentially 
to the bulk value as we move away from the edge. 



(-1)' 



/(r)-(r^ + 3) 
/(r) + (r2+3) 



(107) 



These expressions imply that the correlation lengths (^string) 
for the string correlation are different from ^spin for the spin- 
spin correlation function. 



The infinite-distance limit of the string correlation is given 



as: 



stnng 



(,.2 + 1) (^2 + 9) 



(108) 



It is easy to check that when ?- = eg. (1108b reproduces the 
value 4/9 of the spin-1 AKLT model^. The results are plot- 
ted in Fig.fTTltogether with the correlation length ^string (f)- In 
contrast to the first case |sVBS-I) (see Fig. |9]l, the r oo 
limit of OsJJ^ija is zero since spins disappear from the state 
IsVBS-II) in this limit. 
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FIG. 17. The infinite-distance limit of the string coiTelation function 
(S'l exp[i7r X]J=r^^ 'S'|]S|4.,i) (n oo) as a function of r. The 
value of string coiTelation smoothly decreases from the AKLT value 
4/9 to (no spins left). 



VII. SYMMETRY-PROTECTED TOPOLOGICAL ORDER 

Though the string-order parameter captures the diluted 
Neel order of the Haldane phase, the string-order itself is frag- 
ile under small perturbation&SiS?. Recently, Li and Haldane 
proposed^ to use the structure of the low-lying part of the 
entanglement spectrum (the logarithm of the eigenvalues of 
the reduced density matrix for either of the two partitioned 
systems) as the signature of topological order inherent in the 
state. Pollmann et alji^^ii have investigated the relation be- 
tween the level structure (e.g. degeneracy) of the entangle- 
ment spectrum and discrete symmetries of the system; they 
showed that, for odd-5 spin chains, the existence of (at least 
one of) the three discrete symmetries (time-reversal symme- 
try, link-inversion, and Z2 x Z2 symmetry) guarantees (at least 
two-fold) degeneracy in each entanglement level, while for 
even-S* spin chains, the existence of the above discrete sym- 
metries tells nothing about degeneracy. By this observation, 
they have argued that the Haldane phase in odd- 5 spin chains 
is a stable topological phase protected by discrete symmetries. 

Such arguments can also be applicable to the stability dis- 
cussion of the Haldane-like phase of the present SUSY spin 
models. For instance, the type I sVBS states contains the 
UOSp(l |2) superspin-M multiplet that consists of two SU(2) 
spin multiplets whose spins differ by 1/2. By partitioning 
a superspin-7\/ sVBS infinite chain to two semi-infinite seg- 
ments, there appear two SU(2) spins A//2 and (Af — l)/2 
on the "edge" of each of two sVBS chain segments (hence 
(2A/ + 1) edge states instead of (5 + 1) ones in the usual 
spin-S* VBS states). It is noted that, regardless of the par- 
ity of the bulk superspin M, the sVBS state accommodates 
a half-integer SU(2) spin on the edge. Therefore, for any 
integer-superspin sVBS states, the entanglement spectrum al- 
ways contains a sector consisting of at least doubly degenerate 
levels which come from the half-integer SU(2) spin sector of 
the entanglement Hilbert space. For example, the entangle- 
ment spectrum of the M = 2 sVBS state consists of a dou- 
bly degenerate level corresponding to the doubly degenerate 



fermionic sector and a bosonic level with three-fold degener- 
acy. In fact, we can show that if one of the discrete symmetries 
(link inversion and time-reversal) is present in the SUSY spin 
chains, there is always a sector in the entanglement spectrum 
each of whose levels is at least doubly degenerate. This im- 
plies that the 'Haldane phase' is stabilized regardless of the 
parity of the bulk (integer) superspins. We will report the de- 
tails elsewhere. 



VIII. SUMMARY 

In the present paper, we have constructed a supersymmetric 
extension of the matrix-product states (sMPS) for two differ- 
ent types (I and II) of supersymmetric VBS (sVBS) states and 
exactly evaluated various physical quantities. The sMPS con- 
structed here contains the fermionic elements as well as the 
usual bosonic (i.e. commuting) ones and this slightly compli- 
cates the treatment (for instance, instead of the trace, the su- 
pertrace is used for the periodic systems). We investigated the 
hole-doping behaviors of various correlation functions (spin- 
spin and superconducting) and the spin- and the hole excita- 
tion spectrum. 

In the charge sector, the type I sVBS chains exhibit insulat- 
ing behavior at zero and infinite concentrations of the doped 
holes and the superconducting order parameter is finite only 
for finite doping. In the spin sector, the type I sVBS chains 
interpolate between the usual VBS state and the inhomoge- 
neous VBS state (in the simplest case, it reduces to the MG 
dimer state) at the two extremal limits of hole-doping r = 
and r = 00, respectively. The single-mode approximation 
has been applied to obtain the spin- and the charge excitation 
spectrum. There are two types of low-lying excitations, i.e. 
the triplon and the spinon-hole pair, created respectively by 
the bosonic and fermionic generators of the super Lie algebra. 
The spinon-hole pair is peculiar to the sVBS states; it simulta- 
neously possesses the property of the spin- 1/2 spinon and the 
unpaired hole in the superconducting background. We have 
found that the spinon-hole pair can be the lowest excitation in 
some parameter region of the hole-doping. 

As another class of sVBS states based on a larger (Af^2) 
SUSY, we have introduced the type II sVBS states. In the 
high-doping limit (r 00), the superspin- 1 {M ~ 1) type II 
sVBS state reduces to the totally uncorrected hoIe-VBS state, 
while it reproduces the spin VBS state in the zero-doping 
limit. The type II sVBS state displays qualitatively similar 
behaviors in the spin- and the charge properties except that 
now physical quantities are more symmetric with respect to 
the point r = 1 reflecting that the model contains the equal 
numbers of bosons and fermions. 

We have demonstrated the existence of a hidden order in 
the sVBS states (both type I and II) by calculating the non- 
local string correlations. What is remarkable is that the string 
correlation revives upon hole doping although it vanishes in 
the pure-spin limit r — )- when the spin S = M is even 
integer This may be understood as an example of symmetry- 
protected topological order in SUSY spin chains. 

Though the present work is restricted to ID chains, the 
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sVBS states themselves can be formulated on any lattice in 
arbitrary dimensions, and may generally exhibit resonating- 
valence-bond (RVB) features at finite hole doping. For in- 
stance, an il/ = 2 sVBS state with three species of holes sim- 
ulates the Rokhsar-Kivelson RVB^ in the high-doping limit. 
Such higher dimensional analyses are interesting both theo- 
retically and experimentally, and may be carried by a super- 
symmetric extension of the tensor network method. 
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Appendix A: A crash course on supersymmetry 

1. UOSp(l 1 2) and UOSp(2 1 2) 

The superalgebra UOSP(l |2) consists of the following five 
generators 



S"^ = {ah + b'^a)/2 , Sy = {a)h - h'^a)/{2i), 
S"" = {a)a-h^h)/2 (bosonic) 

K, = ^{x-^fa^ + xpb) 
K2 '^{^^'^ f^'^ ~ ^Po) (fermionic) 
satisfying the (anti)commutation relations: 

[S"' , 5''] = ie^^'S*" (a,6,c = x,2/,z) 



(Ala) 



(Alb) 



(Ale) 



At this stage, the parameter .t, which defines a one-parameter 
deformation of UOSp(l |2), is arbitrary. The second equation 
implies that the fermionic generators K\ and K2 span a two- 
dimensional spinor representation of SU(2). 

Any irreducible representation of UOSp( 1 1 2) is specified by 
superspin l{~ 0, 1/2, 1, . . .). A convenient way of construct- 
ing a superspin-Z representation is to use Schwinger operators 



(bosons a, b and fermion /) satisfying 

na+nb + nf = ah + 6^5 + /V = 2Z (e Z) . (A2) 
Then, the Casimir operator C is calculated as: 

C = S^ + {K1K2 - K2K1) = l{l + 1/2) . (A3) 
The SU(2) subalgebra depends only on a and h: 
= {{na + Tibf + 2{na + nb)] S{S+l) 



S ^ {na + nt)/2 



(A4) 



Since ri/ = 0, 1, a (4?+l)-dimensional superspin-Z represen- 
tation splits into two SU(2) irreducible representations: 



(1)5" = / (n/ 0) • • • (2Z + l)-dim 
(ii) S = 1-1/2 [nj = !)■■■ 2l-Aim , 



(A5) 



which are connected to each other by the fermionic generators 
A'l 2. For instance, the five states in the I = 1 representation 
are: 

(i) 1+) = Ivac), |0) = ajfellvac), |-) = i^f |vac), 



(ii)|t>=al/;|vac), |;)=6j//|vac) 



t fh 



(A6) 



In constructing the sVBS states, we identify (ii) as a one-hole 
state. The ^ = 1/2 case is relevant in realizing the so-called 
superqubi^. 

A two-site system can be treated in the same manner as in 
SU(2); we just define S'°* = S^^' + S'^^ Kf}^ = K^^l + K^^ 
and the corresponding Casimir operator by 

= C(i) +C(2) + 2 {s(i).S(2) + (A7) 
= C(i)+C(2)+25W-5(2) . 

The Clebsch-Gordan decomposition is simply given as: 

1®1~Q®]^®1®---®{21~1/2)®21. (A8) 

So far, the deformation parameter x is arbitrary. However, in 
order for {a\b\ — fojoj — Tf\f\) to behave as a UOSp(l|2)- 
singlet, = r is required. 

By flipping the relative signs of the first and second terms 
in A',, dAlbl ). one may define "new" fermionic operators: 



Di = ]^{-x~^fa^ +xfh) 
D2 = -^ix-'fb^ +xpa). 



(A9) 



The type I sVBS states are not invariant under the transfor- 
mation generated by D^. (Thus, the largest symmetry of the 
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type I sVBS states is UOSp(l|2).) With inclusion of D^, the 
UOSp(l|2) generators satisfy the UOSp(2|2) algebra 



[Sa, Sb] = i^abcSc , {K^.K^} = - (eaa)/ii/S'a , 



[5*0, A'p] = -((Ta)^pi^,y , [5*0, I?^] = -{(Ja)iyfiD„ , 



(AlO) 



where T is defined by 

r = a^a + 6^& + 2/t/ • (All) 

Appendix B: A quick recipe for matrix-formalism 

In this section, we extend the standard formalism for 
bosonic matrix-product states so that we can handle fermionic 
states as well. 

1. Norm 



We begin with the computation of the norm of |MPS). 
Since we consider cases where Aj is made up with both 
bosonic- and fermionic states, a special care has to be taken 
and we proceed step by step. If we write the matrix indicies 
expUcitly, |MPS) reads: 



|MPS)(„,^) 

= ^Ai(a,/3i)A2(/3i,/32)... 



(Bl) 



where the arrow indicates how the order of matrix multiplica- 
tion and the site indices (1, 2, . . . , are related. If the parent 
Hamiltonian H — j+i is defined in such a way that 

hj,j+i{Aj (E) Aj+i) ~ (for all matrix elements) , (B2) 

the matrix indices are physically related to some zero-energy 
degrees of freedom localized at the boundaries {edge states). 
It is important to keep the order (^•) of the string of matri- 
ces. If we adopt the following convention for the hermitian 
conjugation of fermionic operators: 



then the dual of |MPS) reads 

(MPS|(„,^) 

= ^A2(;3l_i,7)..- 

{ft} 

= ^4(7Jl-i)... 

{ft} 

' (B4) 

where A* and denotes a matrix obtained by replacing 
|-) ('I in Aj and its transposition. 

For a periodic chain, the fermion sign has to be treated care- 
fully. Using the identity ipj^TZtpi = STi{TZipiip'i) (the super- 
trace STr is defined in such a way that extra minus signs are 
multiplied for the fermionic sectors. See (|42b) and (|97t , for 
instance), we can express the supersymmetric MPS (sMPS) as 

L-l 

|sMPS)pBc = l[{^pl^z^^+l)^t{'^^lr\o) 



Grassmanii even 



STr|(7^V'l) n(V'!^V'^+i 



(B5) 



= STTiAiA2---AL). 



(/1/2 • ■ • fl-lfi)^ = fifi-1 ■ ■ ■ /2/1 



(B3) 



Since the overlap A*{f3j^i, /3j)Aj{/3j^i, /3j) is a commut- 
ing c-number {transfer matrix), it is straightforward to show, 
by proceeding term by term from the inner most overlap to the 
outer, the following equation: 

(MPS IMPS) = {Ti • • • • • • TL}ic.,a-nn) > (B6a) 

where 

Ti(a,a;^i,/3i) = (a, (a, 
Tl(^l-i,/3l-i;7,7) = A2(;3l_i,7Ml(/3l-i,7) 

rj(/3,_i,/3j„i;^j,/3,) = ^*(/3j-i,/3jM,(/?,-i./3j) • 

(B6b) 

For the purpose of calculating various correlation functions, it 
is convenient to consider generalized overlaps of the following 
form: 

(„,^)(MPS|MPS)(^,5) = {Ti • • - T, • • •Ti}(,,^^^,5) 

- j(a,7;/3,5) , 

which are not necessarily proportional to 5a.-ySp,s for finite- 
L. If the periodic boundary condition is imposed, the norm 
corresponding to the bosonic MPS eq.( l32bb reads 

(MPSIMPS)PBC = J2 {^'^}(c.,/3;o,« = TrT^ . (B8) 

a, 13 

In the case of sMPS, the above expression should be replaced 
with eq.( l46ab . 
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2. Correlation functions 



It is straightforward to obtain: 



Having established the way of evaluating overlaps, it is 
straightforward to extend it to correlation functions. For sim- 
plicity, we only consider bosonic operators here (we will gen- 
eralize the calculation to fermionic operators as well). 

Let us consider first the ordinary two-point correlation 
function: 



_ (MPS lO^Of IMPS) 



(MPS IMPS) 



(B9) 



Since the two physical operators and Oy are bosonic, the 
calculation goes in almost the same manner as in the case of 
norms except that here we have two new matrices: 



Tf0y_l,Py_,-Jy,(3y) = A; {^y _ ^ Jy) O ^ Ay {(Sy _ ^ , ^ y) 

(BIO) 



instead of Tx and Ty. Then, by using 
eq.( |B9l ) may be expressed as: 



the numerator of 



f jnn.-lj,0^j,y-x-lj.O«yL-i/| ^ (Blla) 

I- J (Q,a;7,7) 



Therefore, the matrix-product expression of the correlation 
function is given by: 



(CxCf )(q,7) 



(a,Q;7,7) 



^ J (a., a: 



(Q,a;7,7) 



(Bllb) 



In physical applications, we will encounter the following 
string-like correlation functions: 



n ^ 

yj=X+l 



Oi 



(a, 7) 



(MPS|0;^ of) 0,^|MPS ) 

(MPS IMPS) 



7) 



(B12a) 



Oi 




(a, 7) 



J (a.a; 



(Q,a;7,7) 



(a,a;7,7) 



(B12b) 



In order to calculate the so-called string correlation function 
(see section UVb . we should take: 



, o' 



O^ 



(B13) 



When we consider the expectation values involving 
fermionic operators, the calculation is slightly more compli- 
cated as we have seen in section llllDI 

Appendix C: MPS for Type I sVBS with general M 

The construction of MPS for the 7\f = 1 type-I VBS state 
in section HII Al can be readily generalized to general A/. To 
this end, it is helpful to note that 

(CD 

Since each term on RHS can be written in terms of matri- 
ces with dimensions A/ + 1 or A/, the valence-bond operator 
on LHS may be expressed by a block-diagonal (2A/ + 1)- 
dimensional matrix (a generalization of TZ in (|35]|). Following 
the same steps as in eqs.(l36b and ( |37] |. we obtain the follow- 
ings: 



AMj) - T^{a],b]J])T^{a],b],f])\y), , 



(C2) 



where the polynomials {a — 1, . . . , 2Af+l) are given 
by 

' f 1 Net — 1 / F! _,Q!— l„,J\f— (q— 1) 

for a = 1, . . . , A/ + 1 
fora Af + 2,.. 



, 2A/ + 1 
(C3a) 



and 



for/3 = 1,...,A/ + 1 



^r^M-iCp^iM+2)X^^''+'~>-^yf'-^''+^^z 

for ^ = A/ + 2, . . . , 2A/ + 1 

(C3b) 



with the standard binomial coefficient 
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